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A general spectral transformation is proposed and described"! 

Its spectrum can be interpreted as a Fouriei spectrum or a Laplace 
spectrum. ll;nce it describes the relationship between these two 
transformations. Using this transformation, general time func- 
tions, e.g. functions with exponential growth in both negative and 
positive time directions, can be represented. This general 
spectral transformation is characterized by the introduction 
of two complex frequency variables (p and q) for the positive 
and the negative time domain, resnect ive 1 y . Thus it is possible 
to express the spectral distributions of the Fourier transforma- 
tion (the Dirac function and its derivatives) by rational func- 
tions. The laws and operations of the new method are discussed 
in comparison with the known transformations. Kith the general 
spectral transformation it is possible to define the stability 
criteria for general (i.e. causal and non-causal) systems. As 
an example, its application is shown for the solution of linear 
differential equations, taking into account the initial state 
due to past system excitation. 

1 ^ ^t rodiict ion 

Spectral transformations have gradually become the most 
important calculation procedure in information technology. Unfor- 
tunately, there are two procedures, the Fourier transformation 
(1) and the Laplace transformation [2], each valid under different 
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preconditions. It is also impossible to limit oneself to only 
one of these methods. Of course, the Laplace transformation is 
always applicable when causal signals arc involved (time func- 
tions only for t > 0) and when one wants to describe realizable 
systems. It then has the advantage of yielding the spectrum by 
an algebraic expression of the complex frequency p (in many 
cases a rationally broken function!. 


With the increased use of the system theory introduced by 
K. KUpfmllllcr (3|, however, it is also necessary to consider 
non-cuusal signals and non- rca 1 i zable (idealized) systems. The 
Laplace transformation cannot be used for this, and one must 
resort to the Fourier transformation. When working with statis- 
tical signal theory, one must also used the Fourier transforma- 
tion. Of course, for stationary or semistationary processes the 
Fourier spectrum contains distributions (i.e. the Pirac function 
and its derivatives), which is an undesirable trait that often 
impedes mathematical handling. There also arises the question 
of whether there are different spectra for the same time signal, 
and what these differences may be. This question yields the 
further question of whether there might be a still more general 
spectral transformation that would include both the Fourier and 
the Laplace transformation, and from which both processes would 
develop as special cases. Such a transformation must on the one 
hand not be limited to causal signals, and must on the otiicr 
hand permit representation of all the processes accessible to 
the Laplace transformation, e.g. functions that develop exponen- 
tially. 

/3b4 


The following describes such a general spectral transforma- 
tion and discusses its laws and peculiarities. In particular, 
it also establishes the connection between the Fourier and 
Laplace transformations and permits conversion of their spectra. 
The present report can contain only a short sketch of the proce- 
dure. The reasoning behind the general spectral transformation 
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given in detail in (4); the necessary proofs and the derivation 
of their laws are also contained there^. 

2. Common Grounds an d Di f fere nces b e t ween the Fo u H e r and L aplace 
Transformations 


Both procedures are based on the l-ourier integral: 


4 «• 


fr(/) . 

•• 

(«) 

4 M 


(2) 


The differences are essentially rooted in the manner in which 
the convergence of this integral is reached. Namely, with sta- 
tionary (or unilaterally stationary) time procedures, the Fourier 
integral does not converge. 


The Fourier transformation forces convergence by means of 
a "convergence factor", e.g. e ^1^1, yielding 

4 oo 

■"lim J' tt (/)»«- *l'l , (3) 

• -’•o 
+ «o 

• *0 

(In F.q . 4 a convergence factor e was introduced for reasons 

of symmetry.) Thus it is possible to represent stationary pro- 
cesses, but one gets distributions by passing the limit c-*0. 

For instance, for u(t) ■ 1, one gets 




rrU) = *(/), 



where (f) is the Dirac function with the function integral 
*j6(f)df - 1 and «(f) - 0 for f M- 


For the convergence of the Fourier integral the time func- 
tion must be exponentially limited, i.e. it must satisfy the 
condition 

The author wishes to thank his colleague Prof. Bordewijk, Delft 
Technical University, and his assistant W. Wolf, for stimulating 
discussions on this theme. 
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Thus time functions such as u(t) ■ ct^ are 
re represented in the Fourier spectrum by deri- 
Dirac function. 

e transform is limited to causal time processes 
■ 0 for t < 0. By a substitution, p ■ ju ■ j2iTf, 
convergence of the Fourier integral for complex 
n unilaterally stationary processes, Re(p) > 0 
required. The transformation equations are 

» J . JM ' 


transformation here should be chosen large 
integration path running parallel to the imaginary 
n the convergence domain of IJj^(p). In this way 
ly exponential processes can be represented. The 
a (6) can be transformed to a ring integral 
t < 0 and t > 0, using Jordan's lemma, and ran 
ed simply using the residue theorem of function 
suit of evaluation is called a Heaviside develop- 


V(p) uit) u(t) 

A fijcri ** wrh 

iJpIKe -triflsf 



1. Pole position of a single pole in the 
ex frequency plane, and corresponding time 
ion according to Laplace and Fourier trnas 
t ions . 

a. p-plane b. according to 

• c. real 
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To compare the two procedures, we now observe a spectrum 
as in Fig. 1, with a single real pole at point a, i.e. the 
spectrum 

C) 

p — II 

» 

Depending on whether a < 0 , a ■ 0 or a > 0 , the reverse trans- 
formation yields different results via the Laplace and Fourier 
transformations. For a > 0 (Fig. 1, Case 1), the result is 
the same with both methods. For a < 0 (Case 3) it is different, 
because the integration path with the Fourier transformation 
runs along the imaginary axis of the p plane (for real frequen- 
cies) -- i.e. left of the pole -- while with the Laplace trans- 
form it is right of the pole. The exponentially developing 
process appears in the right time domain (t > 0) with the 
Laplace transformation, while it is reverse-poled and in the 
left domain (t < 0) with the Fourier transformation. jhe 
case of a pole directly on the imaginary axis for a ■ 0 (Case 2) 
is especially difficult, for now one must know how the Fourier 
integral is to be evaluated. If one takes Cauchy's main value, 
one gets 2 Sgn(t) with the Four;er transformation. With the 
Laplace transformation one gets the unit step for t > 0. Thus 
one and the same spectrum yields different time functions with /365 
the two procedures, or in other words, the "spectrum" as defined 
up to now is ambiguous and one must know which method is being 
used . 

3. Definition of a General Spectral Transformation 

One can now combine both methods and find a common repre- 
sentation, if one assigns the time domains t > 0 and t < 0 
different complex frequency symbols. We choose p for t > 0 

and q for t < 0, and write 


I’AGK 

qUALU’V 
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(H) 
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We get a snectrum > with two differentiated complex fre- 

(p.q) * 


quency variables, p and q. Here the partial spectra 

can be calculated by a right- or left-side Laplace ti 
formation, for 


V{m (I)}. (10) 

(M-v)'- vi{-(-0}. («>) 




•V ’-(q) 
Under these 


For the reverse transformation the integration path J from 

-j® to ♦j« is run so that all 
poles p^ of left of 

it and all poles q of U 
arc right of it 

conditions, p ■ q * x can again 
be established for the reverse 
tr'insformat ion . Fig. 2 shows 
the integration path for U(X,X) 
in the X plane. The p poles 



p^ are shown by the symbol )d 
and the q poles of q arc shown 


Fig. 2. Integration path in 
the X plane with the general 
spectral transformation. 

Key: a. X plane 

with the symbol P< . Since the 
position of the p and q poles is arbitrary, this method can 
represent time processes that develop exponentially on both 
sides. With respect to representabi 1 ity , the general spectral 
transformation thus combines the possibilities of the Fourier 
and the Laplace transformations. 
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Fig. 3. General spectral transformation; time function and 
spectrum of a single pole at the frequency zero point. 

Fig. 3 shows the time function, general spectrum and 
Fourier spectrum for a single pole in the frequency zero point, 
and Fig. 4 shows the same for a multiple pole. One can see 
that the distributions occurring with the Fourier spectrum 
(as well as their derivatives with the multiple pole) can be 
represented by rationally broken functions in p and q. 

One gets the Fourier spectrum from the general spectrum by 
the limit passage 

bV(/) ■= lim f f. V — f) . (12) , 

« -*0 

For the Laplace spectrum 

I'Ap). ( 13 ) 

applies, i.c. the partial spectrum dependent on q is 
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Fig. 4. General spectral transformation. Time function and 
spectrum of a g-multipie pole in the frequency zero point. 

simply omitted. 

For instance, for the spectrum of u(t) - 1 (direct current), 

«{/) - lim ‘ ‘ ]. (14) 

■ r <l — t I 

/366 

applies. This means that the following correspondence applies 
for the spectral Dirac function (spectral line): 

;» V 

2 

For the Dirac function differentiated p times, likewise 
^ The conversion ‘^^^^(f) ■ (w) applies. 

8 
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Thus spectral distributions of the pth order correspond to 
the difference between p and q poles of the p^l^^ order. 

According to Fig. 4 this includes time functions increasing 
with the power of p. If one only wants the causal portion (for 
t * 0) (Laplace spectrum), one simply omits the q component. 

Through the initial complication of having to use two 
frequency variables p and q, one gains the advantage of always 
being able to distinguish the positi»re and negative time domain, 
even in the spectrum and .'if obtaining the Fourier *^pectrum (if 
it exists) by passing a limit. The advantage of the Laplace 
transformation is maintained, namely spectral representation by 
rational functions in p and q instead of distributions, and thus 
the possibility of applying tlic rules of function theory (residue 
theorem). In addition, one can represent time processes that 
grow exponentially to both sides. Further to be noted is that 
the procedure is significantly more efficient than the already- 
known bilateral Laplace transform, which requires the convergence 
of the spectrum in a strip of the p plane parallel to the ima- 
ginary axis. 


4. The Development Theorem Acc o rding t o Figcn fu nct io ns 


If U, X is given by a rational broke'.i function, one can 
get a development theorem correspond to the Heaviside theorem 
of the Laplace transformation. First one must split off an 
entire function in p or q, setting p ■ q ■ X. The remaining 
function component is developed into a partia. fraction 

series with the coefficients A 
order u 


IV 


of the p poles p of the 
vg V 

One proceeds correspondingly with the function compo- 


nent gets a partial fraction development with the 


coefficients B of the q poles q 
vg ^ 


of the order u 
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Here is the time derivative of the temporal Dirac 

function fi(t). The function 


Y(t) 


I 


1 for t > 0 
0 for t < 0 


is the unit step. Aside from the distributions in the time zero 
point one gets a development according to the eigenfunctions at 
the po’e points, which are exponential oscillations of the polar 
frequencies multiplied by a power of t. The separation of the 
two time domains (p poles for t > 0 and q poles for t < 0) is very 
evident here. 


5. Operations 

Table 1 contains the most important operations of spectral 
transformation, comparing Fournier, Laplace and general spectral 
transformations. For the derivation and proof of the latter 
one must refer to [4]. Some peculiarities should be noted. 

5 .1. Causal and Ac a usal Dif ferentiation an d In te gr'ation 

By distinguishing p from q in the spectrum, not only the 
time domain but also the time direction becomes distinguishable. 
For in the usual limit passage that occurs in forming the temp- 
oral differential quotients, the causal relation (distinction 
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Spectral 
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1 
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differentiation 

J a n ilf 

dp 

U/< dy 

Spectral 

Integration 

i>«) 

(-0 


•'W 

Spectral folding 
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Temporal 
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b. )causal * causal 

c. akausal > acausal 
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between cause and effect) is lost. This has the result that /3b8 
in differential eijuutions from physics time is symmetrical and 
apparently d i i ec t ion 1 ess , which has already led to may contra- 
dictions in understanding the development of events. (In the 
author's opinion it is unnecessary to resort to the Second Law 
of Thermodynamics to >;ive time a direction on the urounds of 
irreversibility.) lor if one forms the differential quotients 
while ohserviiiR the causality requirement, i.e. with the stipu- 
lation that only past function values may be used, then (At > 0) 

U(/) - Mil - ,1/) 


C). 


Iim 

N •)! 


.I< 


(III) 


Considering the displacement theorem 


M(< — for At >0 (20) 

this yields, after the limit has been passed, to the operation 

( 21 ) 


which is designated as a causal differentiation. Correspondingly, 
an "acausal" differentiation in which the future function values 
are used. 


cri 


> lint 
• If >u 


Mil I 


. 1 /) - Mil) 

.U 


would lead, with the help of the correspondence of the displace 
ment theorem 


M|/ T .h)o " for At > 0 (23) 


to the operation 

• I i\*t \ 


This means that a "causal" d i f ferent i a t io.n requires the multi- 
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plication of the spectrum by p, and an "acausal" differentia- 
tion requires multiplication by q. Thus a t'»mporally causal 
differential equation yields a spectral power expression only 
in p, and thus the conversion factor S of a causal system is 
also only a function of p: 

(W) 

The impulse response s(t)o — • thus a causal temporal 

function. By requiring a causal differentiation, the direction 
of time established, which would otherwise not be the case. 

The equivalent situation applies to integration, in which 
past values of u(t) are used in the causal case and future values 
in the acausal case. Accordingly, the spectrum must be multiplied 
by i (causal) or ~ (acausal). 

5.2. Multiplication and Fo l ding 

In multiplying two general spectra one should note that 
there are no restrictions on the products of p-partial spectra 
and q-partial spectra. However, with a mixed product of p and 
q partial spectra one must require a parallel to the imaginary 
axis, so that all p poles of the one factor arc left of the axis 
and all q poles of the other are right of it. The folding of 
the corresponding time functions corresponds to the product of 
the spectra. 

On the other hand, the folding of the spectra corresponds 
to the products of the time functions. From this it follows 
that folding a p-spectrum with a q-spectrum must always yield 
zero : 

Spectral folding may be defined by a complex integral as follows: 


13 



Vlip.q)* (’tip.q) (*") 

- . ‘ ^ /A (i. A) A. y- A) .1 A. 

Here the integration path Jj of Uj is used, and one should choose 
Re(p) > Re(X) ♦ Re(p 2 ^) and Re(q) < Re(X) ♦ Re(q 2 ^) in performing 
the integration. 

Using spectral folding one can calculate partial spectra, 
if they are not already available: 

. (2H) 

r.(7)- ‘)- (2«) 

This corresponds to a multiplication of the time function by 
Y(t) or Y(-t). 

5.5. Limit Value Theo rems 

If u(t) at t ■ 0 has no distribution, one can calculate 
the left and right-hand initial values for t ■ +0 or t » -0 
from the partial spectrum by the following limit passages: 

«(+(»)=« liiii I /»• r’,(/i)|, CU)) 

HHt>) -»om 

m(~ 0) = lim [ — . (31) 

a* 


Likewise, for the asymptotic final values at t » ♦» or t ■ 
(if these exist) , 

m( , oo) n liin [;t (\ (/<)] , (il-) 

, *. «(— oo) -« lim(— 
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5.4. Hilbert Transformation 


The Hilbert transformation of the Fourier spectrum 


can be expressed especially simply using the general spectrum. 
Here 

u{t)o^ LHq), ( 3 *) 

M(l) o 0,U) & V, ,p) - . (3.^) 


/569 

applies. One thus gets the Hi Ibert -transformed Up(f) from Up(f) 
by reversing the sign of the q partial spectrum and multiplying 
by j. The correspondence sign ^ indicates that the Fourier 
spectrum is obtained from the general spectrum by limit passage 
(Eq. 12). 


6 . Stability Criteria for General Syste ms 

The general spectral transformation permits statements on 
stability for general -- i.e. non-causal -- system functions 
such as occur in system theory. Since the relationship S(p,q) 
•-0 s(t) applies, where S(p,q) is the system function (conver- 
sion factor) and s(t) is the impulse response, the following 
three stability criteria apply: 


a) If s(t) is restricted by energy, i.e. 

oo. 

“O* • 

then all p poles of S(p,q) must be in the open left half of the 
frequency plane and all q poles must be in the open right half. 
The imaginary axis (including oi « ®) is polelcss. 
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b) If s(t) is restricted by power, i.e. 

I * Tit 

r / l'(0|*«W<oo. 

, r .« / ' 

then the p poles of S(p,q) must be in the closed left half of 
the frequency plane and the q poles must be in the closed right 
half, i.e. including the imaginary axis. On the imaginary axis 
(including u» ■ ®) the poles are single. 

c) If s(t) is exponentially restricted, i.e. 

V 1*(0| «l/< oo 

*• m 

for all positive t , then all p poles of S(p,q) must be in the 
closed left half of the frequency plane and all q poles must be 
in the closed right half. Multiple poles may occur on the 
imaginary axis (including u ■ ®) . 

(Note: Such systems develop, for instance, in series 

circuits of ideal integration amplifiers.) 

In addition to system theory, which uses idealized but 
not immed„ately realizable acausal systems (for which s(t) ■ 0 
does not apply for t < 0) , such systems also occur practically, 
if one replaces time t with place x. This is the case in 
"homogenous layer system theory"(5), which serves to describe 
neuron systems (nerve networks in a layer-like structure). Here 
the same stability criteria apply with respect to the local 
stimulus transmittance of such systems. The same problems arise, 
for instance, in realizing digital recursive filters for 
picture processing. 

7. Solution of Differential Equation s with Initia 1 Values 

Here the Laplace transformation is actually particularly 
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well suited, since it can take account of the charge condition 
of the system as expressed in the initial values. On the other 
hand it cannot include charging that occurred in the past. Here 
the general spectral transformation establishes the relationship. 



Fig. 5 . Time graphs of the compulsive force 
Uj(t) and the effect U2(t). - charging 

of the system by u, (t) and homogenous solution 
u,j^(t) for t > 0 ; * ■ compulsive force 

Uj^(t) and its effect U22(^) for t > 0 . 

In Fig. 5 the dotted line represents the charging process 
of the system and the solid line represents the compulsive force 
Uj^(t) and its effect U2^(t) ■ 02jj(t) U22(f)> The term 02|j(t) 

is the homogenous solution of the differential equation for t > 0 
taking into account the intial state. U22(^) is the effect of 
the compulsive force Uj_^(t). Both components should be calcu- 
lated separately. From the differential equation of the system 
(which is assumed to be causal), 

• 

( 3 * 1 ) 

r* 0 r 0 


one gets by transformation into the spectral domain using the 
general spectral transformation: 

» " **. 

, 1 «r 2 (/».</)”> b, p» i\ {p,q). (37 ) 

»-U tm\) 


This yields 




Here the causal system function 
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m j n 

- i. ■ (30) 

r II ' *-« , 

is determined by the (real) coefficients of the presumably 
causal differential equation (36). With the partial spectra 
of Uj , corresponding to both time domains, 

• • * ’ 

• * 1- ('/) . (^'') ■ 


one gets 

I'iip.i) " •'''(/') • f'ltip) i- •\(/') • t'»-('.') . (■•>) 

The first term in Hq . (41) describe the effect of the compul- 
sive force Uj^(t) commencing at t ■ o, so that 

utr.{i)o • ♦(/•). (4:J) 


The second term of Fq . (41) describes the influence of the 
past, i.e. the effect of the cause Uj_(t). Of this, only the /'370 
t > 0 component is of interest -- namely the homogenous solution 
of the differential equation u,^j(t). We must thus calculate the 
p-partial spectrum from the spectral product us ing iE"q . ( 28) , 
and get 

«2ii(/)o •l.v(/i) • ri-(f/)| * . (4:{) 


This spectral expression represents a causal system function H(p) 
describing the homogenous solution (because of the charging): 

//(/») ^ |.s’(/.) r,_(v)i* ' . (44) 

/' 

Thus the replacement circuit in Fig. 6 applies, in which the 
effect of the compulsive power starting at t * 0 and the homo- 
genous solution are represented separately. The latter is pro- 
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duced by a Dirac function 4(t) over H(p) acting at t ■ 0. 


bq.(44) can also be 
transformed into the time 
domain, yielding for the 
homogenous solution 

•‘111(0 - (01 -y(/). 

Here an ambiguity is permis- 
sible insofar as totally 
different functions Uj (t) 
lead to the same function 
U 2 ||(t). Kith the determina- 
s established unambiguously. 
H(p) desciibes the charge condition of the system at the time 
t ■ 0. On the other hand, H(p) can also be calculated using 
the Laplace transformation. If one applies the differentiation 
theorem of the Laplace transformation, which takes into account 
the initial condition, one gets, for t > o, Uj^(t) ■ 0 as the 
compulsive force disappears (cf.(4]). 

■ I ! • 

( 4 ( 1 ) 

i »**0 * r 0 

« 

Here the denominator polynomial is the same as that of S(p) and 
the coefficients of the numerator polynomial are given by the 
initial values of the homogenous solution for t ■ ♦O and its 
derivatives. 

o- ( «). (•‘7) 

f-u 



Fig. 6. Replacement circuit 
diagram of the differential 
equation for the effect of the 
compulsive force U 22 (t) (for 

t > 0) and the homogenous solution 
u,|^(t) based on the charging of 

the system. 


tion of H(p), however, the procc*«s i 


Since U 2 j^o — •H(p) applies, these initial values may also be 
determined in reverse from ll(p) via the limit theorem (30): 
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If one knows the coefficient c^, whether by knowing H(p) from 
bq.(44) in development towards rq.(4b), or from the initial 
value from l'q.(470, one can get the analog computer realization 
for the solution of the differential equation as in Fig. 7. 
(Hero the real constant k is an arbitrary time factor and t 
is the time constant of the integration amplifier used, deter- 
mined by U 2 ■ (i/Ujdt.) 



Fig. 7. Realization of the differential 
equation in analog computer technology with 
integration amplifiers of the system function 
l/pt and a time scale factor k. The coeffi- 
cients a, and b are given by Fq.{36) o’* (39), 
and the Coefficients c^ are given by ’q.(4b) or (47). 

The initial state of the system, which is expressed in the 
coefficients c and which must be known for a realization by 

V 

analog computer, can thus be determined from the system's past 
using the general spectral transformation according to Fq.(44) 
and ( 46) . 
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